Centrality Dependence of Charm Production from Single Electrons Measurement in 

Au+Au Collisions at ^/sNN = 200 GeV 
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The PHENIX experiment has measured mid-rapidity transverse momentum spectra (0.4 < pr < 

4.0 GoV/c) of single electrons as a function of contrality in Au+Au collisions at y^ijvjv = 200 GoV. 
Contributions to the raw spectra from photon conversions and Dalitz decays of light neutral mesons 
arc measured by introducing a thin (1.7% Xq) converter into the PHENIX acceptance and arc 
statistically removed. The subtracted "non-photonic" electron spectra are primarily due to the semi- 
leptonic decays of hadrons containing heavy quarks (charm and bottom). For all centralities, charm 
production cross section is found to scale with the nuclear overlap function, Taa- For minimum-bias 
collisions the charm cross section per binary collision is Nc-/Taa = 622 ± 57 (stat.) ± 160 (sys.) /jb. 

PACS numbers; 25.75.Dw 

In central Au-|-Au collisions at ^sjviv = 200 GeV tt^'s and charged hadrons are strongly suppressed at high 
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transverse momentum (px) In contrast, a mod- 

est high-pT enhancement is observed in d+Au collisions 
at the same energy 0, |^. This strongly suggests that 
the suppression observed in Au+Au collisions is caused 
by final-state effects {e.g., parton energy loss in a dense 
medium produced in the reaction 

Heavy quarks (charm and bottom) are complementary 
probes of the hot and dense matter produced in high 
energy heavy ion collisions. Due to their large masses, 
charm and bottom cross sections are calculable via pQCD 
and their yield is sensitive to the initial gluon density (8| . 
It has been predicted that heavy quarks suffer less en- 
ergy loss than light quarks while traversin g pa rtonic mat- 
ter due to the "dead cone" effect 0, IkI lllj . This can 
be studied through systematic measurements of the pr 
spectra of open heavy flavor. In addition, the open- 
charm yield is an important baseline for understanding 
J/ip production which has been predicted to be either 
suppressed or enhanced 0| in the presence of de- 
confined quarks and gluons. 

The PHENIX experiment observed that inclusive sin- 
gle electrons in central and minimum-bias Au-|-Au col- 
lisions at i/satat = 130 GeV were produced in excess 
of purely "photonic" contributions (primarily due to tt" 
Dalitz decays and conversion of 7r° photons in the de- 
tector material) This excess is consistent with the 
expected charm production, assuming that it scales with 
the number of binary nucleon-nucleon collisions {Ncoii), 
or equivalently, with the nuclear overlap function, Taa- 
In this Letter, we present results on the single electron 
measurement in ^/snn = 200 GeV Au-|-Au collisions. 
The new data have higher statistics and smaller system- 
atic errors than the 130 GeV data, allowing us to measure 
charm production as a function of collision centrality. 

The data used in this analysis were collected by the 
PHENIX experiment during the 2001 run period 
of the Relativistic Heavy Ion Collider. A coincidence 
of the beam-beam counters (BBC) and the zero degree 
calorimeters (ZDC) provides the minimum-bias trigger 
(92.2j:3;o % of the 6.8 ± 0.5 barn Au-hAu inelastic cross 
section). The centrality is determined by the correlation 
between the multiplicity measured by the BBC and the 
energy of spectator neutrons measured by the ZDC. The 
BBC also measures the collision vertex, z, with resolution 
a — 0.7cm. Events are required to have \z\ < 20cm to 
eliminate electrons originating from the central magnet. 

Charged particles are measured by the PHENIX east- 
arm spectrometer (|?7| < 0.35, Acf> — 7r/4) with reso- 
lution cTp/p ~ 0.7% ® 1.0% p(GeV/c). Tracks are re- 
constructed with the drift chamber (DC) and the first 
layer of pad chambers (PCI) and confirmed by requiring 
an electromagnetic calorimeter (EMC) matching hit 
within 2 standard deviations in position. Electron candi- 
dates are required to have at least three associated hits 
in the ring imaging Cerenkov detector (RICH) that pass 
a ring shape cut, and are required to pass a timing cut 



in either the EMC or the time-of-flight detector. After 
these cuts, a clear electron signal is observed as a narrow 
peak at E/p ~ 1. By requiring —2a < (E ~ p)/p < 3cr, 
background from hadrons, which deposit only a fraction 
of their energy in the EMC, and non- vertex electrons, 
which have mis-reconstructed momenta, is further re- 
duced. Remaining background in the electron sample, 
due to accidental coincidences between RICH hits and 
hadron tracks, is estimated (w 10%) and subtracted by 
an event-mixing method. 

Inclusive electrons contain two components: (1) "non- 
photonic" - primarily semi-leptonic decays of mesons 
containing heavy (charm and bottom) quarks, and 
(2) "photonic" - Dalitz decays of light neutral mesons 
(tt'^, t], ?/, p, Lu and (p) and photon conversions in the 
detector material. To separate these two components, 
a photon converter (a thin brass tube of 1.7% radiation 
length surrounding the beam pipe at r = 29 cm) was 
installed. 
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FIG. 1: Shown vs pt (a) Raw spectra measured with 
the converter in (open circles) and out (closed circles), (b) 
Ratio of the converter in/out yields {Rcn, points) and 
ratio of photonic yield with/without the converter {R^, 
line and shaded band), (c) Ratio of non-photonic to photonic 
yields (Rnp, points) and contribution from kaon decays 
(dashed line). 



We analyzed 2.2 M (2.5 M) events with the converter 
in (out). The corresponding raw electron pT spectra for 
minimum-bias collisions are shown in Fig.^a). The pho- 
ton converter multiplies the photonic contribution to the 
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electron yield by a factor R^: 



n: 



Conv — out 



n: 



Conv — in 



R^N2 +(l-e)-iV: 



(1) 
(2) 



Here ivf""^"™ (^Cont,-o«t^ ^-^^ measured electron 
yield with (without) the converter; N2 (iV^°"-T) is the 
electron yield due to the photonic (non-photonic) com- 
ponent; and e(« 2.1%) represents a small loss of elec- 
trons due to the converter. We next define Rcn as 
the ratio of the raw electron yield with and without the 
converter. Dividing Eq. (O by Eq. ^ and defining 
Rnp = N2°'^-'^/N2, one has: 



Rcn 



n9' 



onv — tn 



(1 - e)R 



NP 



j^Conv-out 1 + RnP 

If there were no contribution from non-photonic compo- 
nent {Rnp = 0), then Rcn — R-y- 

The photonic electron yield per photon is approxi- 
mately given by Y (x where S is the Dalitz branch- 
ing ratio per 7 relative to 2j (for 71^,77,77') or I7 (for p, 
uj and (f)) decay, and t is the thickness of the conversion 
material in radiation length {Xq). The factor | is the ap- 
proximate probability for a 7 to convert in one Xq. Plug- 
ging in S""" = 0.6%, t « 1.1% {t « 2.8%) for converter in 

(out) we find R^ = yCon«-»rYyCor,i,-o«t _ ;^_g rj.^^^^ 

is some pt dependence in the complete formula for Y 
and the value of 5 is species-dependent {6^ w 0.8%), so 
we perform a full GEANT simulation with and with- 
out the converter to calculate R^. We determine R^ for 
7r° and rj separately. We use the 7r° spectrum measured 
by PHENIX as the input for the tt" simulation and as- 
sume rriT scaling (px — > y^p^T-KM^^^^M^, normalized at 

high Pt to ri/ir'-' ~ 0.45 ±0.1, based on the world data of 
tj/tt^ ratio) to obtain the input for the ij simulation. Con- 
tributions from other mesons which undergo Dalitz decay 
{rj' , p,u!,(l)) are small (6% at px = 3GeV/c, and smaller 
at lower px)- Since they have S we assign them 

R^ = When calculating the combined R^ we use the 
particle ratios at high px [r]' /ir^ = 0.25 ± 0.13, /o/tt^ = 
w/tt" = 1 ± 0.5, (j)/iT° = 0.4 ± 0.2). The (j)/TT° ratio used 
here is consistent with our tt*^ and (f) measurement |l8j . 
The uncertainties in the particle ratios are included in 
the systematic uncertainties of R^. For this method it is 
essential that the amount of material is accurately mod- 
eled in the simulation. We compared the yield of identi- 
fied photon conversion pairs in the data and in the sim- 
ulation and conclude that the simulation reproduces R^ 
within ±2.0%. This uncertainty is included in the overall 
systematic uncertainty. 

Fig. n^b) shows Rcn and R^ vs. px- Rcn gradu- 
ally decreases with increasing px, while i?^ slightly in- 
creases with Px- The difference between Rcn and R-y in- 
dicates the existence of non- photonic electrons. Fig.Q^c) 
shows Rnp obtained from R^ and Rcn using Eq. ||2Jl. 
Rnp increases with px and is more than 30% for px > 
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0.6 GeV/c. The small amount of conversion material in 
the PHENIX detector allows a sensitive measurement of 
Rcn- 

Background from kaon decays (K — > Trev) and di- 
electron decays of p, lo and </) remain in the non-photonic 
electron yield. The background from kaon decays is es- 
timated with a GEANT simulation using the kaon px 
spectrum measured by PHENIX Jl^ as input. The con- 
tribution of kaon decays to the non-photonic yield, shown 
in Fig. n^c), is 18% at px = 0.4GeV/c and decreases 
rapidly to less than 6% for px > 1 GeV/c. 

To calculate background from the e~^e~ decays of p, uj 
and (j), we first generate spectra by applying mx scaling to 
the PHENIX 7r° spectrum, as described above. The con- 
tribution of these decays to the non-photonic electrons is 
< 3% for all Px- Background from J/i/; e+e" decays 
and from Drell-Yan pairs is negligible. Possible enhance- 
ment of low mass di-leptons through tt + tt — > p — > e+e~ , 
as reported in Pb+Pb collisions at the SPS |2(|], would 
contribute to the non-photonic electrons. However, this 
is neglected since the estimated p contribution in the ab- 
sence of enhancement is only « 0.6% over all px- 

After these backgrounds are subtracted the only other 
significant source of non-photonic electrons is the semi- 
leptonic decay of heavy flavor (charm and bottom). The 
raw spectrum of heavy flavor electrons is corrected for 
geometrical acceptance (e^eo), track reconstruction effi- 
ciency (crec) and electron identification efficiency {eeio) 
determined by GEANT simulation. The efficiency egeo x 
tree is about 11% of dNe/dy, and teiD is about 65% as 
confirmed with electrons identified through photon con- 
version. Correction of multiplicity dependent efficiency 
losses, estimated by embedding simulated electron tracks 
into real events, is pT-independent and increases from 
5% to 26% from peripheral to central collisions. The 
Icr systematic uncertainty of these corrections is 11.8%. 
Fully corrected heavy flavor electron spectra are shown in 
Fig. 121 for minimum- bias collisions and for five centrality 
bins. 

PHENIX has also measured the heavy flavor electron 
spectrum in p-\-p collisions at ^/snn — 200 GeV . The 
lines in Fig. 13 show the best fit curve of this spectrum, 
scaled by Taa for each Au+Au centrahty bin. Here, Taa 
is the nuclear overlap function calculated by a Glauber 
model Q (Tabled). The Au+Au data points are in rea- 
sonable agreement with the p+p fit in all centrality bins. 

To quantify the centrality dependence of heavy flavor 
production, we calculated the integrated yield dNe/dy 
(0.8 < Px < 4.0GeV/c) and fit it to AN^^^, where 
a = 1 is the expectation in the absence of medium ef- 
fects. In this comparison, most of the systematic ef- 
fects will cancel. Figure |21 shows dNe/dy{0.8 < px < 
A.O)/NcoU vs. Ncoii for minimum-bias and five centrality 
bins in Au+Au collisions and p + p collisions. We find 
a = 0.938 ± 0.075(stat.)±0.018(sys.). li p + p data is 
included, a = 0.958 ± 0.035(stat.). This shows that the 
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FIG. 2: Fully corrected heavy flavor electron pr spectra for 
different Au+Au centralities scaled by successive factors of 
ten for clarity. Error bars (brackets) correspond to statistical 
(systematic) uncertainties. Curves are described in the text. 



total yield of heavy flavor electrons for all centralities is 
consistent with Ncou scaling. 
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FIG. 3: Non-photonic electron yield (0.8 < pr < 4.0 GeV/c) 
measured in Au+Au reactions at 200 GeV scaled by NcoU 
(left-hand scale) as a function of centrality given by Ncoii- 
This electron yield translates to the electron cross section per 
A'' A'' collision in the above pT range (right-hand scale). The 
yield in p + p collisions at 200 GeV is also shown IgJ] . 

For each centrality bin we scale the heavy flavor elec- 
tron spectrum [pT > 0.8 GeV/c) by Taa and fit it with a 
PYTHIA calculation of the electron spectrum resulting 
from leading order charm and bottom production. We 
used PYTHIA 6.205 with a modifled set of parameters 



(described in [Ij) and CTEQ5L PDFs [23. Based on 
experimental input [Ull^l we modified the PYTHIA 
default charm ratios, using instead /D^ — 0.45 ± 0.1, 
DJD° = 0.25 ± 0.1, Ac/D° = 0.1 ± 0.05. This gives a 
c-quark e total branching ratio of 9.5 ± 0.4%. The 
scaled charm and bottom cross sections are treated as 
fit parameters, although we find that our data are re- 
stricted to PT values which are only sensitive to charm 
production. We evaluated the systematic error due to 
background subtraction (« 21%) by refitting to the elec- 
tron spectrum at the minimum and maximum of its la 
systematic error band. The change of the pt range for 
fitting the non-photonic electron spectrum gives 3% sys- 
tematic error for minimum-bias collisions. The system- 
atic error due to the PYTHIA spectral shape (« 11%) is 
dominated by the uncertainty in (/c^) — 1.5±0.5GcV/c. 
Different PDFs yield a systematic error of 6.2% for the 
rapidity-integrated cross section. Systematic errors in 
Taa are tabulated in 0. These systematic errors are 
added in quadrature to give the overall systematic error 
on the charm cross section. For minimum-bias collisions 
we obtain xII^^^^\v=o — 143 ± 13(stat.) ± 36(sys.) /ib 
and N^s/Taa = 622 ± 57(stat.) ± 160(sys.)^b. Results 
for all centrality bins are shown in Table U The STAR 
collaboration reports somewhat larger charm cross sec- 
tion (1.3 ± 0.2 ± 0.4 mb per NN collision) in p + p and 
d+ Au colhsion at y^s^N = 200 GeV Hi] . The next-to- 
leading order pQCD calculation of charm cross section is 
300 to 450 ^lh j2Q]. 



TABLE I: Centrality bin, number of NN collisions, nu- 
clear overlap function, charm cross section per NN collision, 
and total charm multiplicity per NN collision, in -y/sjviv = 
200 GeV Au+Au reactions. 



Cen- 


N^oii 


Taa 


1 d^c- 1 

Taa dy \y=o 


N^-c/Taa 


trality 




(mb-i) 




i^^h) 


min. bias 


258+25 


6.14+0.45 


143+13+36 


622+57 +160 


0-10 % 


955+94 


22.8+1.6 


137+21+35 


597+93 +156 


10-20 % 


603+59 


14.4+1.0 


137+26+35 


596+115+158 


20-40 % 


297+31 


7.07+0.58 


168+27+45 


731+117+199 


40-60 % 


91+12 


2.16+0.26 


193+47+52 


841+205+232 


60-92 % 


14.5+4.0 


0.35+0.10 


116+87+43 


504+378+190 



It should be noted that final-state effects only influence 
the momentum distribution of charm; they have little or 
no effect on the total open charm yield. Therefore, our 
results indicate Ncoii scaling of the initial charm produc- 
tion, as expected for point-like pQCD processes. pQCD 
calculations without charm quark energy loss and hydro- 
dynamic calculations assuming complete thermalization 
of charm quarks predict very similar heavy flavor elec- 
tron spectra for pT < 2GeV/c jl^. Differentiating be- 
tween these opposite physical pictures is only possible for 
Pt > 2.5GeV/c, where statistics of the current analysis 
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are limited. 

In conclusion, we have measured single electrons from 
heavy flavor decays in An + Au collisions at y^s]v7V = 
200 GeV. We observe that the centrality dependence of 
charm quark production is consistent with Ncoii scaling, 
as expected for hard processes. The much larger Au + 
Au data set collected by PHENIX in the 2003-04 run will 
allow us more detailed exploration of medium effects on 
heavy quark production, both through deviations of the 
heavy flavor electron spectrum from Ncoii scaling, and 
also through a measurement of charm quark flow. 
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